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Exactly Solvable Model for the Decay of Superdeformed Nuclei 
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The history and importance of superdeformation in nuclei is briefly discussed. A simple 
two-level model is then employed to obtain an elegant expression for the branching ratio 
for the decay via the El process in the normal-deformed band of superdeformed nuclei. 
From this expression, the spreading width for superdeformed decay is found to be 
determined completely by experimentally known quantities. The accuracy of the two- 
level approximation is verified by considering the effects of other normal-deformed states. 
Furthermore, by using a statistical model of the energy levels in the normal-deformed 
well, we can obtain a probabilistic expression for the tunneling matrix element V. 

1. Superdeformed Nuclear Decay 

Superdeformation is one of the most interesting examples of collective phenomena in 
atomic nuclei. Since its original experimental observation in 1986,^ the properties of 
these high-spin rotational bands have fascinated experimentalists and theoreticians 
alike. When combined with precisely measured branching ratios and decay rates, 
a thorough theoretical understanding of the mechanism by which superdeformed 
(SD) nuclei are formed and then decay into normal-deformed (ND) bands promises 
to provide a window into nuclear structure unlike any other. 

One of the major barriers to using SD decay to understand nuclear structure 
has been that there is no clear way to link the quantities measured in experiment 
directly with those which might shed light on the internal dynamics of the nu- 
cleus. Despite enormous strides made by experimentalists to measure observables 
precisely,^''^'^ theorists have failed to reach a consensus on just what to do with these 
data. Ideally, we require a model which, while accounting for the rich physics of the 
SD nucleus, is also as simple as possible to allow for easy extraction of quantities 
of interest to theory. 

The typical SD nuclear experiment^'^'^'^ creates nuclei in a high angular mo- 
mentum state in the SD potential well. These nuclei then lose rotational energy 
by E2 transitions, eventually reaching a low enough angular momentum that it 
becomes energetically favorable to decay to states of the same angular momentum 
in the ND well. As each SD nucleus continues to decay, more strength moves into 
the energetically preferred ND band rather than continuing down the band of SD 
states. In practice, most of the decay-out of the SD band happens over the space 
of only two or three SD states, after which essentially all the strength has moved 
into the ND band. For a schematic diagram of this process, see Figure 1(a). 
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Fig. 1. (a) Schematic diagram of the decay-out process. The SD band loses essentially all its 
strength to an ND band over the course of a few states, (b) The problem is commonly modeled 
with a double-well potential in deformation space. The shape of the potential is a function of the 
angulcir momentum of the nucleus. 

It is customary to model the decay process as a double well problem, as in Fig- 
ure 1(b). The nucleus is a quantum mechanical system moving in a spin-dependant 
potential in deformation space. The potential consists of two wells, and the eigen- 
states of each well in the absence of the other make up the pure SD and ND bands. 
In such models, it is the shape of the potential at various angular momenta which 
contains information about the underlying nuclear structure. Thus, an important 
goal for theories of this type is that they provide a method to relate the potential 
to experimental observables. 

2. Two-State Model 

The two-state model for SD nuclear decay ^ is an exactly solvable approximation to 
the decay-out scheme outlined in the previous section. The key assumption is that 
only one ND level participates significantly in the decay of a particular SD level of 

the same angular momentum. This reduces the problem to the familiar two-level 
mixing scenario of introductory quantum mechanics. The great advantage of this 
approach is that the two-state model rigorously yields simple formulae connecting 
experimental observations to the physics of the potential barrier. 

Keeping only the ND level of the same angular momentum nearest in energy 
to the decaying SD state, one can simplify the system to that shown in Figure 2. 
The dynamics of the model are thus completely determined by four parameters: 
the difference in the unmixed levels' energies A = ejv — £s, the energy widths for 
coupling to the electromagnetic environment Fg and F^v, and the tunneling matrix 
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Fig. 2. Schematic diagram of the two-level model for SD nuclear decay. Es and £jv are the energies 
of the unmixed SD and ND states, respectively; and Tff are energy widths due to interaction 
with the electromagnetic field; and V is the matrix element for tunneling through the barrier. 

element V. A fifth parameter, the mean level spacing in the ND band Dn, can be 
taken as the unit of energy. 

In the absence of the tunneling parameter V, the system reduces to two isolated 
Breit-Wigner resonances, and the retarded Green's function in the energy domain 
is simply 



1 







Go{E) = ; I (1) 

\ ^ E-A+irN/2 I 

in the 5, N basis. Coupling between the two levels is given by the simple perturba- 
tion matrix 

The full Green's function is given exactly to all orders in V by Dyson's Equation. 
Its inverse is 



(3) 



The retarded Green's function contains all information about the time evolution 
of the system. We are specifically interested in the branching ratios, which are given 
by Parseval's Theorem to be 



where % is either S or A''. The ND branching ratio follows: 

rjvrV(riv + ri) 



(4) 



(5) 



Vs + rjvrV(riv + ri) ' 

a result we first noted in Ref. 7. Here F^^ is the spreading width for tunneling 
through the barrier: 

r* = (6) 

A2 + r 
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Tabic 1. Values of and (V) for various SD decays in the A Ri 150 and A fa 190 
mass regions. (F^) from Rcf. 2 is tabulated for comparison. Values for _FjVi Tg, 
Fjv, and Djv are also from Rcf. 2. is calculated from Eq. (7), either directly 
or, in the case of '^^^Dy(26), statistically, as discussed in the text, while {V) is 
given by Eq. (15). Angular momenta for the decaying SD states are given in 
parentheses. 
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152Dy(28) 


0.40 


10.0 


17 


220 


11 


35 


41 000 




0.81 


7.0 


17 


194 


78 


87 


220 000 


i94Hg(12) 


0.40 


0.108 


21 


344 


0.072 


5.0 


560 


"4Hg(10) 


0.97 


0.046 


20 


493 


1.6 


35 


37 000 



where T = ^"+^^ . 

In this case, Equation (6) is the exact result of Eq. (4), and it is also what one 
expects from correct appHcation of Fermi's Golden Rule.'' Wc emphasize that, in 
general, ^ tli6 ensemble average of F-'-. In much previous literature (e.g., 

Ref. 6), it was assumed that the two are equivalent, giving erroneous results. In 
fact, they are equal only in the continuum Umit of overlapping resonances.^ This is 
clearly not an appropriate approximation to most nuclei of interest, in which Djv 
exceeds Fjv by several orders of magnitude. Indeed, in real nuclei the discrepancy 
between F^- and its ensemble average can be as much as three to four orders of 
magnitude (see Table 1). 

Equation (5) demonstrates the power of using a simple, intuitively understand- 
able model. We can instantly see that the result is correct, since it can be read as 
the branching ratio of a two-stage decay problem. That is, is the probability for 
the nucleus to tunnel through the barrier and then to decay down the ND band. 
Furthermore, it is eminently useful: we see that, given a set of experimental results, 

is known. Inverting Eq. (5) we find 

, ^ FnTnTs 

Fw-F;v(Fiv + Fs)' 

so that the spreading width is fixed by experimental results. The recent work of 
Wilson and Davidson^ has shown how to relate F^- (Eq. (7)) to the barrier height 
by assuming it is equivalent to a fusion-like tunneling rate. 

Since F^- is a positive definite energy width, Eq. (7) implies a lower bound in 
the two-level model for the decay width due to El transitions in the ND band, 

FnTs 

This prediction of the two-level model is especially significant due to the relatively 
high uncertainty in current values for Fjy. 

The estimated value^ of Tn for ^^^Dy(26) violates inequality (8). This could 
indicate a breakdown of the two- level approximation. However, such a conclusion 
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would be precipitous due to the large uncertainty in the current best estimate of 
Fiv, which could be off by a factor of two or more. An estimate of F^^ for this case 
can be obtained in the two-level model by cutting off the normal distribution for 
Fjv with inequality (8). The median value of so estimated is presented in Table 
1, assuming the standard deviation ar^ = Fjv = 17meV. 

3. Gaussian Orthogonal Ensemble 

We have seen that the two-level model gives an expression for the spreading width 

completely in terms of experimentally measurable quantities. Since F^ is clearly 
related to the potential barrier between the two wells, this is a great stride to- 
ward the goal laid out in the first section, namely to find a simple relationship 
between obscrvables and one or more parameters related directly to the shape of 

the potential. 

The spreading or decay width F^, however, is a rate, and less directly related to 
the shape of the potential than V, the tunneling matrix element itself. The problem 
in relating V directly to experiment is clear from Eq. (6): we would need to know 
the value of A, which is unfortunately not accessible with current experimental 
techniques. Lacking a complete solution of the nuclear structure problem, we derive 
a probabilistic theory for V, with its roots in the Gaussian Orthogonal Ensemble 
(GOE). 

The GOE is a "structureless" assumption for the level distribution of the un- 
mixed ND band, in which we assume only that the levels are eigenstates of some 
Hermitian Hamiltonian with time reversal and SO (3) symmetry. The states are 
characterized by the Wigner surmise,^ a probability distribution for the spacing 
sfjv between levels: 

= |se— (9) 

Since we are interested in finding a probability distribution for unmixed levels, 
the SD and ND spectra are uncorrelated. Thus, given s, we have a rectangular 
probability density function for A 

The Heaviside step function specifies that A is the detuning of the nearest neigh- 
bor. 

Application of basic probability theory leads us to the density function for A: 

ViA) = Vs{A)P{s)ds = ^erfc (v^;^) • (H) 

Here erfc(a;) is the complementary error function of x. The mean of this distribution 
is £'iv/4. 
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From Eq. (6) we note that 




(12) 



which imphes 




(13) 



This ahows calculation of the desired probability density function for V: 





(14) 



where A is given as a function of V by Eq. (12). 

We have arrived at a probability density function for V solely in terms of ex- 
perimentally measurable parameters. The mean of this distribution is 



Values of (V) for specific SD decays are given in Table 1. In general, we can expect 
that (V) will be a typical value of V for nuclei measured in the laboratory, since 
the standard deviation of the distribution (14) is comparable to the mean (15). 

Equation (15) is a profound success of the two-level model. In fact, even by 
including more levels, a probabilistic statement like this represents the most infor- 
mation one can have about V with the current types of experiments. 

4. Beyond the Two-level Model 

In the previous two sections, we solved an approximate model exactly. In reality, 
of course, the ND band is semi-infinite, starting at the bandhead and continuing 
upward. We now address the validity of the two-level approximation. 

The simplest way to address this issue is to add a second ND level.'^ The branch- 
ing ratio is now calculated in a three-dimensional Hilbert space, in exactly the same 
way as before. The inverse Green's function is 



whore the subscripts 1 and 2 pertain, respectively, to the nearest and next-nearest 
ND states to the decaying SD with the same angular momentum. 

It is straightforward to obtain an analytic result for the three-level system by 
applying Eq. (4), but since we desire to know how changed the total branching 
ratios are from the two-state results over a range of parameters, it is perhaps more 




(15) 




(16) 
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Fig. 3. Numerical results comparing the two- and three-state total branching ratios. Superscripts 
indicate the model used whereas subscripts represent the nearest (1) and next-nearest (2) neighbor 
in the three-level case, (a) A 190 mass region, as determined by Ts/r^r = 10"'^ (b) A ^ 150 
mass region with Fg = Fjy. In both cases T/D^ = 10~*. The levels were taken at their mean 
detunings within the GOE under the assumption that they lie on either side of the SD state, and 
in both cases Vi =V2 = V. 

useful to calculate them numerically. Figure 3 shows the results for the A « 190 
and A pa 150 mass regions. 

We note that the figure shows remarkably good agreement between the two- and 
three-level total branching ratios in the mass-190 region, particularly for values of 
V that are likely according to Eq. (15). In the mass-150 region, as for ^^"^Dy, they 
may disagree by as much as 25% (7 = 28) to 40% (/ = 26) due to the larger 
value of the coupling matrix element V and the smaller level spacing, i. e., higher 
density of states. This supports the validity of the two-level approximation for the 
mass-190 region, in the sense that the addition of the third level has not affected 
the physically measured quantities too significantly. In the mass-150 region, the 
three-level approximation yields a significant correction to the two-state branching 
ratios for values oIV/Dm ^ 0.15. 

Further support for the two-level approximation is given by Ref. 10, in which 
Dzyublik and Utyuzh exactly solve the problem in the approximation of an infinite 
(in both directions), equally spaced ND band. Their results also show remarkably 
strong agreement with the two- level approximation for typical values of A. 

5. Conclusions 

Our three-state results, together with the results of Ref. 10, demonstrate that the 

two-state model is sufficient to describe the dominant decay-out process of SD 
nuclei. Within the two-state model, we have shown that the decay out of an SD 
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level via the El process in the ND band is a two-step process, whose branching ratio 
(5) is expressed in terms of three measurable rates, Ts, Tn, and Wc have also 
shown how to determine the tunneling matrix clement V (Eqs. (14) and (15)) from 
the measured values of and a statistical model of the ND band. It is hoped that 
these results will help to clarify the nature of the decay-out process in SD nuclei, 
and that their elegance will inspire further promising studies such as Ref. 8. 
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